Scaling predictions for radii of weakly bound triatomic molecules 
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The mean-square radii of the molecules 4 He3, 4 He2— 6 Li, 4 He2— 7 Li and 4 He2— 23 Na are calcu- 
lated using a three-body model with contact interactions. They are obtained from a universal 
scaling function calculated within a renormalized scheme for three particles interacting through 
pairwise Dirac-delta interaction. The root-mean-square distance between two atoms of mass tua 
in a triatomic molecule are estimated to be of de order of C \/h 2 /[m A (E 3 ~ E 2 )], where E 2 is the 
dimer and E$ the trimer binding energies, and C is a constant (varying from ~ 0.6 to ~ 1) that 
depends on the ratio between E 2 and E- 3 . Considering previous estimates for the trimer energies, 
we also predict the sizes of Rubidium and Sodium trimers in atomic traps. 
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PACS numbers: 21.45.+V, 34.10.+X, 36.40.-c 



I. INTRODUCTION 

Weakly bound molecules are large size quantum sys- 
tems in which the atoms have an appreciable probabil- 
ity to be found much beyond the interaction range, and 
at the same time tiny changes in the potential parame- 
ters can produce huge effects in the properties of these 
states pj. The best illustration of such systems is the ex- 
perimentally found 4 He2 dimer 0, with J (r 2 ) = 52 ±4 A 
and binding energy 

£ 2 =1.1 + 0.3/ - 0.2 mK0. Other 
examples of weakly bound molecules are found through 
the experimental realization of Bosc-Einstein condensa- 
tion (BEC) 4], where the possibility to change the ef- 
fective scattering length of the low-energy atom-atom in- 
teraction in the trap to large positive values by using 
an external magnetic field jj, |5j can produce very large 
dimers. In fact, weakly bound molecules in ultracold 
atomic traps were reported in Ref. |fj. The binding en- 
ergy of the 87 Rb2 dimer formed in a Bose-Einstein con- 
densate was measured with unprecedented accuracy 0. 
Ultracold Na2 molecules have also been formed through 
photo-association One should note that in limit of 
an infinite atom-atom scattering length tuned by Fes- 
hbach resonances in the trap, in principle the Efimov 
condition^ can be achieved, in which an infinite number 
of weakly-bound trimers exist. The formation of weakly 
bound trimers in ultracold atomic systems have not been 
reported till now, but recently the recombination coef- 
ficient rate was used to predict trimer binding energies 
of some specific atomic species that are being studied in 
atomic traps [Tof . 

Theoretically it is possible to exist weakly bound 
molecules of zero-angular momentum states in triatomic 
systems, as for example, in the extensively studied He 



trimer system (see, e.g. Ref. 01 an d therein). These 
molecules are special due to the large spatial size which 
spreads out much beyond the potential range [Til IT^| . 
In such trimer, the calculations of the mean-square dis- 
tance of each 4 He atom to the corresponding center-of- 
mass have been performed for the ground and excited 
states fill. [FJ . and also for the mean-square interatomic 
distance [Tj. These sizes are of de order of 5 to 10 A 
for the ground state of the 4 He3 molecule, and of about 
50 to 90 A for the excited Efimov state [Tl|. Therefore, 
the system heals through regions that are well outside 
of the potential range, in which the wavefunction is es- 
sentially a solution of the free Schrodinger equation, and 
where the physical properties of the bound system is de- 
fined by few physical scales. For example, the dimen- 
sionless product of the mean-square interatom distance 
with the separation energ y of one atom from the trimer 
is not far from the unity [ll| in the ground and also in 
the excited state, despite of the large difference between 
such energies E^ /E^p w 50 (E^ is the binding en- 
ergy of the n-th trimer state). So, as already discussed 
in Refs. [l3L IliL Il5l IT^| . we should note that quite natu- 
rally the binding energy is the scale that dominates the 
physics of the trimer. One should remember as well that 
the collapse of the three-body system in the limit of a 
zero-range forced?] makes the three-body energy one of 
the scales of the system, beyond the two-body energv[l3|. 

The calculation of the low-energy properties of the 
three-body system can be performed with a renormal- 
ization scheme applied to three-body equations with 
s— wave zero-range pairwise potential 0, I 11 this 

approach, one can fix the three-body ground-state (the 
three-bo dy p hysical scale), and the two-body scattering 
lengths [13. Consequently, all the detailed information 
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about the short-range force, beyond the low-energy two- 
body observables, are retained in only one three-body 
physical information in the limit of zero-range interac- 
tion. 

In the present work, we first study the mean-square 
distances of one atom to the center-of-mass (CM) sys- 
tem and between two atoms in the ground and excited 
states of triatomic molecules of the type 4 He2— X, where 
X= 4 He, 6 Li, 7 Li and 23 Na. Next, using trimer energies 
derived from the recombination coefficient rates we 
make estimates of the corresponding sizes of Rubidium 
and Sodium trimers. We introduce and calculate scaling 
functions that describe the different radii as functions 
of the physical scales of the triatomic system obtained 
in the limit of a zero-range interaction. In this way, we 
are generalizing the concept of scaling function, that was 
previously introduced in Refs. [l^ and to study the 
behavior of bound and excited virtual Efimov states Q 
in terms of triatomic physical scales. 

The scaling function depends only on dimcnsionless 
ratios of the binding energies of two and three atoms, 
and the ratio of masses of the different atoms. In that 
sense our conclusions apply equally well to any other low- 
energy triatomic system. The validity condition for the 
scaling relations is that the interaction range must be 
small compared to particle distances, which is the case 
for weakly bound three-body systems. 

The paper is organized as follows. In section II, we 
present the Faddeev equations for the spectator functions 
for a triatomic system with two equal particles a and a 
third one /?, and the form factors from which the different 
mean-square radii are obtained. Also in this section we 
discuss the gen eralization of the scaling function defined 
in Refs. (T5L Il8l | to describe the different radii. In section 
III, we present our numerical results for the mean-square 
distances of one atom in respect to the CM system and 
between two atoms in the ground and excited states of 
triatomic molecules. Our conclusions are summarized in 
section IV. 



II. RENORMALIZED THREE-BODY MODEL 
AND FORM FACTORS 

In this section, we introduce the generalization of the 
scaling function defined in Refs. and [1^, to De used 
to obtain the different radii. We write down the cou- 
pled renormalized equations for the spectator functions 
and the expressions for the form factors which allow the 
calculation of the different mean-square distances. 



A. Subtracted Faddeev Equations 

Throughout this paper we use units such that h = 
m a = 1. For a = 4 He, ^ 2 / m4 Hc = 12-12 K A 2 . After 
partial wave projection, the s— wave coupled subtracted 
integral equations, for two identical particles a and a 



third one /?, are given by 
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The mass number A is given by the ratio m^jm a . The 
plus and minus signs in @ and |0J refer to virtual and 
bound two-body subsystems, respectively. 

In the present context that we have three particle sys- 
tems with two identical ones, it is worthwhile to call 
the attention to two particular definitions of three-body 
quantum halo states: the Borromean states [T^. where 
all the two-body subsystems are virtual (a — a — /?); and 
the tango states |2(| , where the a — (3 subsystems is vir- 
tual and aa is bound (aa — (3). Note that the virtual 
pair of particles is denoted with a dash between the sym- 
bols. The Borromean case corresponds to positive signs 
in front of the square-root energy of the subsystems in 
both Eqs. © and J3J, implying in the weakest attrac- 
tive kernel of Eqs.QJ an d © among all the possibilities 
of signs in the two-body scattering amplitude. And, for 
the tango three-body system, we have negative sign only 
in front of ^e aa in Eq.Q, with positive sign in front of 
^Je a p in Eq. (J3J| . So, a more effective attraction occurs in 
a tango state than in a Borromean case. Of course that, 
if all the two-body subsystems are bound, the effective 
attraction is maximized; and, if all such subsystems are 
unbound (virtual), the effective attraction is minimized. 

One can extend the classification scheme of three-body 
quantum halo states of the type aa/3 |2l|. considering 
the four possibilities, for increasing values of the mag- 
nitude of the effective attraction in Eqs. and J3J. 
The weakest attractive situation corresponds to the pre- 
vious defined Borromean-type (only virtual subsystems) 
(a — a — P). The tango situation (aa — 0) is followed by 
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an three-body system with a — a virtual and a/3 bound, 
that we represent by (a /3a) halo system. Three-body 
system with the strongest effective attraction has all the 
subsystems bound and it is represented by (aa/3). 

We solve Eqs. Ill Kill in units such that the three-body 
subtraction point /i( 3 ) is equal to one|l8|. The corre- 
sponding dimensionless quantities are: €3 = E$/iOh\, 
e aa = £ m /^ 3) , e Q/3 = E a/3 /iJ, 2 ( 3 y The three-body physi- 
cal quantities can be written in terms of the three-body 
binding energy E 3 when first the value of fi3 3 *. is deter- 
mined from the known value of S3. Therefore, the results 
for the renormalized model appear when the subtraction 
point energy is written as a function of £"3 and conse- 
quently the three-body quantities naturally scale with 
E 3 . Finally, the scaling functions are obtained when the 
dimensionless product of physical quantities are written 
as a function of the ratios between two-body energies and 
E3. 



B. Scaling functions for the radii 

The existence of a three-body scale implies in the low 
energy universality found in three-body s yste ms, or cor- 
relations between three-body observables [la. l22j . In the 
scaling limit [bH ITsj ■ one has 



O (E, E3, E aa ,E a p, ) (£3 



A VE/E 3 , ^/E aa /E 3 , , jE aP /E 3 , A , (7) 



where O is a general observable of the three-body system 
at energy E, with dimension of energy to the power r\. In 
the present paper we discuss only the situation that we 
have only bound subsystems (aa/3); however, the anal- 
ysis could be easily extended to other three-body halo 
systems, as the Borromean, tango and (a /3a) systems. 

In the case of the mean-square separation distances, 
(r 2 ) with 7 = a or f3, i.e, the distance of the atom 7 to 
the CM; and (r„ 7 ), i.e, the distance between the atoms 
a and 7, the scaling functions are of the form: 



\J (r 2 )S 3 = TZ 7 (^Je aa /e 3 , , \Je a p/e 3 , A^j 



and 



\j ( r l 7 ) 5 '3 = K al (^Je aa /e 3 , , ^Je a(3 /e 3 , A 



(8) 



(9) 



where £3 is the smallest separation energy of the three- 
body system, i.e., S3 = min (E 3 — E aa , E 3 — E a p). Two 
particular situations are worthwhile mentioning, one is 
the case of trimer systems (A = 1), where the scaling 
functions above reduce to: 



S 3 = TL 



(V e 2/e 3 
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and 



(11) 



The other special situation is found for e Q7 = where the 
dimensionless product of the square radii and triatomic 
binding energy depend only on the mass ratio: 



^)E 3 = TZ 7 (A) , 



and 



y/ir&Es = K ai (A) 
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C. Form factors 

The mean-square radii are calculated from the deriva- 
tive of the Fourier transform of the respective matter den- 
sity in respect to the square of the momentum transfer. 
The Fourier transform of the one and two-body densities 
define the respective form factors, Fp(q 2 ) and F ai (q 2 ), 
as a function of the dimensionless momentum transfer q. 
For the mean-square radius of the particle 7 (= a or j3) 
to CM, we have 
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Fp{q 2 ) 
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And, for the mean-square distance between the particles 
a and 7, we have 
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where 



F a ~ ( (q 2 ) = j d 3 yd 3 zy al (y,z + l)y a , ( (y,z~l) . (17) 

The above triatomic wave- functions in momentum space 
are given in terms of the spectator functions Xa-f- 



^ aa (y, z) 



X [x«a(M)+Xa3(l^-|l)+ Jfa/j(|*+||) 



^<*(3{y,z) 



(18) 



£3 + 



A+X J2 



z 1 + 



A+2 tfl 



1 + 



A+1^2 



A+2 $2 



• ( xaa(W- j?rrl) + Xo t gi\y\)+XcxB(\z+ y 



A- 



4 



where z is the relative momentum of the pair and y is the 
relative momentum of the spectator particle to the pair 
in units of /i(3) = 1. Note that the sub- indices of ^ in 
Ea. (|18|l just denote the pair of Jacobi relative momenta 
used to evaluate the wave-function. For cry with 7 = a 
or /3, one has the relative momentum between a and 7 
and the relative momentum of the third particle to the 
center-of-mass of the system cry. 



valid for the ground state, and consequently deviations 
in the scaling plot are stronger for this state. 
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III. RESULTS FOR TRIATOMIC RADII 

Our analysis has considered some particular three- 
body molecular systems, in which the three-body ground- 
state energy and the corresponding energies of the two- 
body subsystem is known theoretically for 4 He trimer nil, 
4 He 2 - 6 Li, 4 He 2 - 7 Li and 4 He 2 - 23 Na [23|. In Ref. |nj, 
the authors have considered realistic two-body interac- 
tions; their results for the ground and excited states radii 
are appropriate for our purpose of comparing with the 
present scaling approach. 

The ground and excited Efimov state energies of the 
4 He3 molecule were extensively studied in the scaling ap- 
proach of Refs. 0, with results that are in very 
good agreement with realistic calculations. This lead us 
to conclude that other details (beyond the dimer and 
trimer ground-state energies) presented in the realistic 
interactions, that have been used, are quite irrelevant 
to the existence of Efimov states. These features val- 
idates a universal scaling function, relating the trimer 
ground-state, the dimer and the weakly bound excited 
three-body energy state. Realistic calculations for the 
excited states of He trimer approaches reasonably well 
the scaling limit 0,0], which suggests to investigate the 
scaling limit of other observables like the different radii 
defined in Eqs. (|14|l and 1)16(1 . The conditions for the va- 
lidity of the present approach are that the atoms should 
have a very shallow and short-ranged two-body interac- 
tion and the binding energy close to zero, i.e., the ratio 
between the interaction range and dimer size should be 
much smaller than 1. These are indeed the cases we are 
considering. 

The results for the radii of He3 molecule in the ground 
and excited state are shown in Fig. 1, in th e form of a Q&iyjfi'* /(E^rm-Rc), in the limit of zero pairwise binding 
sc aling plo t. The dimensio nless pr oducts y/(r%)S~3 and energies. Our results for (r 2 a(s ) are 0A5/E 3 for N = and 
V ( r aa)^3 as functions of y/E 2 /E 3 are shown in the fig- 
ure and c omp ared to the realistic calculations, obtained 
from Ref.jlJ. Our calculations for the ground and ex- 
cited state are practically the same, which would be the 
case if the energies in respect to /1 2 are in fact going to 
zero, i.e., the scaling limit. The results for \J (r^}S 3 and 
)S 3 for the excited state are in good agreement 
However, for the ground state 
The excited 



FIG. 1: The dimensionless products y (ra)S3 (lower curves) 

and \J {raa)Ss (upper curves) as functions of yj E2 / E$ ■ Our 
results for the ground state and first excited state are shown 
respectively, by solid and dashed lines. Realistic calcula- 
tions from Ref.jllf . for yj (r^fsl are given by empty squares 
(ground state) and empty circles (excited state); and, for 

\J (yaaiSi, by full squares (ground state) and full circles (ex- 
cited state). 

In Fig. 2, the results for the dimensionless products 
V( r la) E 3, \J(rlf3}E 3 , and yJ(r 2 )E 3 , as functions of A = 
m/3/mon for E aa = E a p = are shown. We perform 
calculations for the ground (N — 0) and excited (N = 1) 
states, as indicated in the figure. One observe in the 
upper frame of Fig. 2, that the results almost saturates 
above A w 3 to the values found in the limit of A = 00. 
The calculations for A — 00 give for (r 2 Q ) the values of 
0.69/.B3 for N=0 and 0.61/.B3 for N = 1. Therefore, for 
the ground states, the root mean-square distance between 
two 4 He in the triatomic molecules can be estimated by 



0.40/i?3 for N — 1. The saturation value for y/(r^ a )E 3 is 



with the realistic result, 
the results show a deviation of about 20%. 
state size is about ten times larger than the correspond- 
ing size of the ground state. Therefore, the scaling limit 
is better approached in the excited state, which is much 
larger than the interaction range, which is not strictly 



achieved fast with increasing A than for y ( r ap) E 3 which 

depends on the difference in the masses of the atomic 
pair. The mean-square distance of one of the atoms 7 
(= a or (3) to the center-of-mass of the molecule can 
be obtained from the lower frame of Fig. 2, where we 



plot J (r^)E 3 as a function of A. One sees that, for the 

infinitely heavy /3-atom, the results for (r 2 ) are the same 

of the (r^p), while yj \r^)E 3 — as the heavy particle 

should rest in the CM of the molecule in this limit. We 
remind the reader that the ratio between the binding 
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TABLE I: Results for different radii of the molecular sys- 
tems aafil, where a = 4 He and f3 is identified in the first col- 
umn. The ground-state energies of the triatomic molecules 
and the corresponding energies of the diatomic subsystems, 
obtained from Ref. |23| . are given in the second, third and 
forth columns. (r 2 7 ) is the corresponding mean-square dis- 
tance between the particles a and 7 (= a,/3). (r~) is the 
mean-square distance of 7 from the trimer center-of-mass. 






E (0) 
^3 




E a p 












(mK) 


(mK) 


(mK) 


(A) 


(A) 


(A) 


(A) 


4 He 


106.0 


1.31 


1.31 


9.45 


9.45 


5.55 


5.55 


6 Li 


31.4 


1.31 


0.12 


16.91 


16.38 


10.50 


8.14 


7 Li 


45.7 


1.31 


2.16 


14.94 


13.88 


9.34 


6.31 


23 Na 


103.1 


1.31 


28.98 


11.66 


9.54 


8.12 


1.94 



FIG. 2: For the triatomic aa/3 system, with 7 = a, (3, it is 
shown the dimensionless products y {r 2 7 )i?3 [upper (a) plots] 
and <J (ry)E3 [lower (b) plots], as functions of A = mpjm a , 
in the limit E aa = E a p — 0. \J (r 2 ) is the root-mean-square 

distance of particle 7 from the center-of-mass, and y (r 2 7 ) 
is the root-mean-square distance between the particles a and 
7. The results for the ground-state (N = 0) are shown with 
solid line (7 = a) and dot-dashed line (7 = /3); and, for the 
excited state (N — 1), with dashed line (7 = a) and dotted 
line (7 = P). 



energies for N = and N = 1 is about 500 for A = 1 [£j 
and 300 for A — 10, while the dimensionless products of 
square radius and energy changes only around 10%. 

In the above, we have considered examples of molecules 
with two helium atoms. However, our results presented 
in Figs. 1 and 2 are more general, such that we can ex- 
tend the estimates to other atomic systems. Of particular 
interest is the analysis of possible formation of molecu- 
lar systems in experiments with ultracold trapped gases. 
By considering, for example, the estimates of trimer en- 
ergies obtained from the recombination coefficient, given 
in Ref. [Tpj , within our approach we can predict the cor- 
responding trimer sizes. In this case, we have a = 
and A = 1 in the previous equations and in the figures. 
However, our unit for a specific trimer of an atom with 
A nucleons will be h 2 /m A = (48.48/A) K A 2 . 

In Table I, we present our results for the different 
radii of the ground state (N — 0) of the weakly bound 
molecular systems aa(3, where a = 4 He and j3 = 4 He, 
6 Li, 7 Li and 23 Na, obtained from the known theoretical 
values of E^°\ E aa and £a/3[I3- Our calculation for 
4 He 3 of y/(r^ a ) gives 9.45 A which is 14% off the value 
11 A obtained in the realistic variational calculations of 
Ref. The same quality of agreement is found for 

(r 2 ) which in our calculation is 5.55A compared to 
6.4A of Ref. The quality of the reproduction of the 
realistic results by our calculations are quite surprising 



TABLE II: Results for the size of trimer systems predicted 
in Ref. 01- The sizes are given by (r 2 ), the root-mean- 
square distance between the atom a and the center-of-mass 
of the trimer system. The atoms of the trimer are identified in 
the first column. For each dimer energy, given in the second 
column, we have two possible trimer energies (columns 3 and 
5), with the corresponding radii given in the columns 4 and 
6. The trimer estimates, given in Ref. [l(|, for 87 Rb|l,-l), 
are for noncondensed (*) and condensed (^) trapped atoms. 



Atom 


E 2 


E 3 




E'3 






(mK) 


(mK) 


(A) 


(mK) 


(A) 


23 Na|l,-l) 


2.85 


7.75 


12 


3.06 


38 


87 Rb|l,-l)* 


0.17 


0.56 


22 


0.175 


114 


87 Rb|l,-l)t 


0.17 


0.47 


25 


0.183 


91 


85 Rb|2,-2) 


1.3 x 10" 4 


2.4 x 10" 4 


1293 


1.7 x 10" 4 


1944 



in view of the simplicity of the present approach, where 
the only physical inputs are the values of the dimer and 
trimer binding energies. The several different radii of the 
molecules 4 He2- 6 Li, 4 He2- 7 Li and 4 He2- 23 Na have values 
larger than those found in the 4 He3 which makes plausi- 
ble that our predictions are even better in quality. 

We point out that the results in Table I, for the 4 He 
dimer sizes inside the molecules shrink in respect to the 
free value of 52A, due to the large values of the trimer 
binding energies. Qualitatively this is explained just by 
considering that the dimer size scales roughly with the 
inverse of the square-root of its binding energy inside 
the molecule, which can be estimated to be 2/3 of the 
molecule binding, from which one finds for that the dimer 
has sizes around 10 A, close to the ones we have found in 
Table I. 

In Table II, we are also presenting results for differ- 
ent radii of the trimers predicted in Ref. |To| . from 
where we obtain the energy of the dimer and the most 
weakly bound trimer energies of 23 Na,\F = l,mj? = —1), 
87 Rb|F = l,m F = -1), and 85 Rb|F = 2,m F = -2), 
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where \F,mp) is the respective hyperfine states of the to- 
tal spin F. We are presenting the mean-square distance 
from each atom to the center-of-mass of the correspond- 
ing trimer. From Fig. 1, one can also obtain the corre- 
sponding mean-square distance between the atoms. We 
observe that one value of the recombination rate is con- 
sistent with two values of the most weakly bound trimer 
energy, as discussed in Ref. 10]. Therefore, we present 
two possible values for the radii that are consistent with 
the corresponding weakly bound trimer energies. Ac- 
tually, we should also mention that, in a trap, one can 
achieve dimer and trimer molecules with very large sizes, 
following the possibility to alter the corresponding two- 
body scattering length 

Finally, it is interesting to recall the results for the 
hyperradius calculations obtained by Jensen and collab- 
orators [2l|. From their scaling plot, one can observe 
that the hyperradius of a tango system is bigger than the 
hyperradius of a Borromean system, for the same three- 
body energy. This point can become very clear, for in- 
stance, if we take as an example the results shown for the 
\H (filled circles in Fig. 2 of Ref . [2ll| 1 . and estimate the 
dimensionless product of observables, (p 2 )mB jh 2 (prod- 
uct of the x-axis and y-axis in Fig. 2 of Ref. |2lj], where p is 
the hyperradius and B the three-body binding energy). 
When going from the tango to the Borromean configura- 
tion, this product will decrease. Therefore, if one keeps 
the same binding energy, the Borromean system would 
be more compact than the tango system. Extending this 
analysis to (a. (3 a) halos, where (a — a) is virtual and 
(a/3) is bound, and also to all-bound pairs (aa/3), one 
should expect that the sizes increase when going from 
Borromean states to halos with all-bound subsystems, 
while keeping the three-body energy fixed. Within our 
approach, the scaling relations are expected to be fol- 
lowed in all the cases, in the limit of a zero-range inter- 
action. The Borromean halo (the less attractive system), 
in order to have the same three-body energy as a tango 
state, should be more compact, in agreement with the 
scaling plot of Ref. [2l|. In a realistic case, our scheme 
is expect to produce better results in the all-bound case, 
when comparing systems with the same three-body en- 
ergies. This occurs because the all-bound case would 
have the most extended wave-function; consequently, the 
range of the potential, in relation to the size, would have 
the smallest value, satisfying better the validity condi- 
tion for the scaling relations, which is that the interaction 



range must be small compared to particle distances. 



IV. CONCLUSIONS 

The mean-square radii of the triatomic molecules, 
4 He 3 , 4 He 2 - 6 Li, 4 He 2 - 7 Li and 4 He 2 - 23 Na are calculated 
using a renormalized three-body model with a pairwise 
Dirac-delta interaction, having as physical inputs only 
the values of the binding energies of the diatomic and tri- 
atomic molecules. Presently, we have considered molec- 
ular three-body systems with bound subsystems, due to 
the available data. The renormalized zero range model 
can also be applied to the cases where at least one of the 
subsystems is virtual. When comparing systems with the 
same three-body energies and potential ranges, he valid- 
ity of the renormalized zero range model is expected to 
be better in the all-bound case, because this case corre- 
sponds to the most extended wave-function. 

The validity of the present framework is substantiated 
by the agreement of our results for the different radii with 
the realistic potential model calculations of Ref.^lJ for 
4 He3 ground and excited states which are within about 
14%. These results are quite surprising in view of the 
simplicity of the approach, where the only physical inputs 
are the values of diatomic and triatomic binding energies. 
We predicted for the first time, as far we know, the values 
of several different radii for 4 He2- 6 Li, 4 He2- 7 Li and 4 He 2 - 
23 Na molecules, from the theoretical values of the binding 
energies calculated in Ref. [23 ■ These other molecules are 
in general larger than the He-trimer indicating that our 
radii predictions for these triatomic ground states can be 
even better in quality than those found for He3. 

In view of the actual relevance of ultracold atomic sys- 
tems that are being experimentally studied, and the pos- 
sibility to observe the formation of molecular systems in 
trapped condensates, we also present results for the sizes 
of rubidium and sodium trimers, considering the binding 
energies that were recently estimated 10] from analysis of 
the corresponding three-body recombination coefficients. 
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